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ABSTRACT 

High-magnetic-field pulsars represent an important class of objects for studying the relationship 
between magnetars and radio pulsars. Here we report on four Chandra observations of the high- 
magnetic-field pulsar J1718-3718 (B = 7.4 x 10 13 G) taken in 2009 as well as on a re-analysis of 2002 
Chandra observations of the region. We also report an improved radio position for this pulsar based on 
ATCA observations. We detect X-ray pulsations at the pulsar's period in the 2009 data, with pulsed 
fraction 52%±13% in the 0.8-2.0 keV band. We find that the X-ray pulse is aligned with the radio 
pulse. The data from 2002 and 2009 show consistent spectra and fluxes; a merged overall spectrum is 
well fit by a blackbody of temperature 186l 18 eV, slightly higher than predicted by standard cooling 
models, however, the best-fit neutron star atmosphere model is consistent with standard cooling. We 
find the bolometric luminosity Lj^ = 4+2 x 10 32 ergs -1 ~ 0.3-E 1 , for a distance of 4.5 kpc. We compile 
measurements of the temperatures of all X-ray detected high-i? pulsars as well as those of low-i? radio 
pulsars and find evidence for the former being on average hotter than the latter. 
Subject headings: pulsars: individual (PSR JI7I8— 3718) — X-rays: stars — stars: neutron 



1. INTRODUCTION 

Over the past few decades, our knowledge about the 
neutron star family has increased significantly. Previ- 
ously, only one kind of isolated neutron star was known, 
the radio pulsars. Powered by their rotational energy, 
they are also called rotation-powered pulsars (RPPs). 
Their luminosities are generally much lower than the ro- 
tation energy loss rate, E. X-ray observations have led 
to the discovery of several new classes of isolated neu- 
tron stars; these include ma gnetars and X -ray isolated 
neutron stars (XINS£3); see IKaspil d20loh for a recent 
review. They all exhibit distinctive properties different 
from those of conventional RPPs. 

Magnetars are isolated, slowly rotating (known periods 
in the range of 2-12 s) X-ray pulsars, having thermal and 
non-thermal X-ray luminosities that are in many cases 
much higher than their spin-down luminosities. Some are 
characterized by repeating X-ray/7-ray bursting activity 
and therefore are called soft gamma repeaters (SGRs). 
Others are less active and are characterized by their per- 
sistent X-ray pulsations; these are classified as anoma- 
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lous X-ray pulsars (AXPs). However, the distinction be 
tween these two classes has been increasing ly blurred, as 
some sources show properties of both ( e.g. iGavriil et al 
[2001 IKaspi et all 12001 IWoods et 
2009bt TKaneko et all 120101: llsrael et alj 
2011). They generally have very high magnetic fields 
(B ~ 10 14 -10 15 G), inferred assuming that their spin- 
down rates are solely a result of magnetic dipole radia- 
tion. It is generally believed that their X-ray luminosities 
are powered by the decay of the ultra- high magnetic fields 
dDuncan fe Thomp son 1992: lThompson fc Duncan]ll995l : 
iThompson fc Duncanlll996t IThompson et al.ll2002f). Fo r 
reviews of m ag netars, see IWoods fe Thompson! (|2006{k 
IKaspil (|2007f ) or IMereghettil (j2008H 

XINSs are a small group of slowly rotating (known 
periods in the range of 3— 11 s), nearby (dis t ance < 
500 p c) neutr on stars (see IKaspi et ali 120061 IHaberj 
120071 iTurollal 120091 for reviews). Emitting appar- 
ently thermal X-ray spectra, they show no hard X- 
ray emission. No radio counterparts have been found 
for these neutron stars. Given their long periods and 
expected small beaming fractions, it is possible that 
their radio beams are m isaligned with our line of sight 
(|Kondratiev et all 12009ft . Therefore, it is not clear 
whether or not they are intrinsically radio quiet. Tim- 
ing observations of XINSs have revealed relatively high 
inferred magnetic fields (~l-3 xl0 13 G), and spin-down 
ages o f the order of 10 6 years (IKaplan fc van Kerkwiikl 
20051: IZane et all l200l Ivan Kerkwiik fc Kaplan! 120081: 
Kaplan fc van Kerkwiikl [2009al) for some of them. Un- 
like magnetars, they show no bursting activity and are 
much less luminous. However, their X-ray luminosi- 
ties are comparable with their spin-down power and sig- 
nificantly higher than those of normal RPPs of sim- 
ilar ages (IKaplan fc van Kerkwiikl I2009a|) . Therefore, 
IKaplan fc van Kerkwiik (2009a) suggest that the cool- 
ing of XINSs is likely affected by magnetic fiel d deca y 
heating as predicted in theory by I Arras et"all (|2004l) . 
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iPons et all (|2007l ) and lAguilera et al.l (j2008ft . An alter- 
native explanation is that the XINSs are surroun ded by 
fall-b ack disks and are heated due to accretion (|Alparl 
l200l . 

One likely crucial group of pulsars for understanding 
the relationships between RPPs, magnetars and XINSs 
is the high-magnetic-field RPPs. There are now sev- 
eral known RPPs that have spin-down magnetic fields 
close to or higher than those of magnetars. Some of 
them are radio pulsars. Sharing properties with both 
classes, these high-B pulsars could be transition objects 
between RPPs and magnetars. Indeed, some magne- 
tars are now known to emit at radio wavelengths, and 
magnetar-like bursting behavior has been seen in one 
high-B pulsar. XTE J 1810-197is a transie nt AXP, first 
detected in outburst (jlbrahim et al.1 120041 ) . This mag- 
netar, originally not emitting in the radio band, was 
observed to h ave radio pulsation s one year after its X- 
ray outburst (|Camilo et al.l l20(36h . Also, the magnetar 
IE 1 547.0—5408 shows radio pulsations ()Camilo et al.l 
2007). Though not a radio pulsar, the high-P» rotation- 
powered (B = 4.9 x 10 13 G) X-ray PSR J1846-0258 
exhibited a sudden, magnetar-like X-ray outburst that 
lasted for a few weeks in 2006 dGavriil et al.l 120081 : 
IKumar fc Safi-Harbl I20M INg et al.l 120081 ) . Thus, it is 
possible that the high-B RPPs are magnetars in qui- 
escence. Recently, a new magnetar, PSR J1622— 4950, 
was discovered via its active radio e mission, yet is rel- 
atively X-ray-faint (|Levin et al.ll2010D . Another magne- 
tar, SGR 0418+5729, detected via its bursting activities, 
was found to have magnetic field B < 8 x 10 12 G, well 
below that of the other magnetars, suggesting that a 
strong surface dipole magnetic fiel d might not be nec- 
essary for magnetar-l i ke behavior (Ivan der Horst et al.l 
120101: iReaet al.l 120101: lEsposito et al.ll2010D . These dis- 
coveries further suggest that there could be a large, un- 
seen population of quiescent magnetars, some of which 
may be 'disguised' as radio pulsars. 

PSR J1718— 3718 is a rad io pulsar discovere d in the 
Parkes Multi-beam Survey (jHobbs et al.ll2004D . It has 
period P = 3.3 s and spin-down rate P — 1.5 x 10 -12 . 
These imply a characteristic age r c = P/(2P) = 34 
kyr, spin-down power E = 4tt 2 IP/P 3 — 1.6 x 10 33 
ergs~ x where / = 10 45 gcm 2 is a fiducial moment of 
inertia of the pulsar, and surface dipole magnetic field 
B = 3.2 x 10 19 (PP) 1 / 2 G = 7.4 x 10 13 G, which is 
the second highest of all known RPPs and is higher 
than that of AXP IE 2259+586E3(B = 5.9 x 10 13 G). 
PSR J17 18-3718has a disp ersion measure (DM) of 373 
cm^pc (jHobbs et all 12004ft. Based on t he DM and 
the NE2001 model ((CordesjLLaziol [200l . the best- 
estimated distance to the pulsar is ~4.5kpc. However, 
the NE2001 model provides a poor estimate of pul- 
sars' distances when they are near the Galactic center 
(jGaensler et al.1 120041 ) . Indeed the distances estimated 
based on NE2001 for pulsars in the nearby clusters NGC 
6221 and NGC 6403 are a factor of ~2-3 smaller than 
their true distances (|Gaensler et al.l l2008f) . Therefore, 
we suggest that the true distance of PSR J1718— 3718 is 
probably in the range ~4.5-10 kpc. 

An X-ray source was serendipitously detected at the 

10 http : //www .physics .mcgill . ca/$\sim$pulsar /magnetar /main 



radio position of PSR J1718-3718 in a 2002 Chan- 
dra observation. iKaspi fc McLaughlin! (|2005l ) found that 
this X-ray source had a soft, thermal-like spectrum, and 
therefore is the likely X-ray counterpart of the radio 
pulsar. However, due to the limited photon statistics 
(see Table [1]), the spectral results were not very con- 
straining. Also the coarse time resolution (3.24 s) in the 
timed exposure mode observation prevented any pulsa- 
tions from being detected. Deeper Chandra observations 
with higher time resolution were proposed and conducted 
in 2009. Interestingly, a large period glitch occurred 
between 2007 September and 2009 January (Manch- 
ester & Hobbs, 2011 in preparation). Four ChandraX- 
ray observations, each separated by ~ 2 months, were 
taken in the hope of detecting X-ray variability, possi- 
bly associated with the glitch, as occurred in the 2006 
outburst of PSR J1846 -0258 (|Kuiper fc Hermsenll2o5l 
iLivingstone et al.ll2010L e.g.). 

Here we report on a temporal analysis of the four new 
Chandra observations of PSR J1718— 3718, as well as on 
a spectral analysis which also includes the archival 2002 
observation. 

2. OBSERVATIONS AND RESULTS 

Four observations of PSR J1718— 3718 were taken 
with the Chandra X-ray Observatory in 2009. Each had 
~33 ks of live time (see Table [1] for details). In these 
observations, the pulsar was positio ned on the Advanced 
CCD Imaging Spectrometer (ACIS. lGarmire et al.l[2003T l 
S3 chip with y-offset of O'.l and Z-offset of 0fl8 from 
the aim point. The other ACIS chips were turned off. 
The data were taken in 1/8 subarray mode (only photon 
events from 1/8 of the CCD were read out in this mode), 
in order to achieve time resolution of 0.44 s, sufficient 
for timing this 3.3-s pulsar. In the 2002 observation, the 
pulsar was detected on the S2 chip of ACIS ~ 8' off the 
aim point, wit h a total of 99 counts in th e 0.8-2.0 keV 
band (Table [Q IKaspi fc McLaughlinll2005l) . 

We started our analysis with the level 2 event files, 
which are the products of the standard reprocessing III 
FT. and analyzed the data using the tools provided in 
CIA£3 version 4.2 (CALDB version 4.2.0). 

2.1. Imaging and Source Position 

PSR J1718— 3718 was detected in all five observations 
using the celldetect tool in CIA0. The best source po- 
sitions as reported by celldetect were slightly differ- 
ent from one observation to another (Table [lj . This is 
the result of the small pointing uncertainty of the Chan- 
dra satellite. We found the average source position to be 
R.A.=17:18:09.83(1) and decl.= -37:18:51.5(2) (J2000), 
where the uncertainties are the standard deviation of the 
detected positions. 

To look for extended emission, we compared the 
image of PSR J1718— 3718 from each observation 
with a simulated point-source image generated by the 
Chandraiay traceirn (ChaRT, a.k.a the Chandra point- 
spread-function[PSF] simulator) and the MAR]£3 tool in 
CIAD 4.2. We used ChaRT to produce a collection of 

11 http: //exc .harvard. edu/ciao/repro_iii .html 

12 http://cxc.harvard.edu/ciao/ 

13 http://cxc.harvard.edu/chart/ 
html 1 14 |http : //space .mit ■ edu/ CXC/MARX/ 
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TABLE 1 

Chandra Observations of PSR J1718-3718 



ObsID 


Date 


MJD 


t f a (s) 


Offset 13 


tf (ks) 


R.A., Decl. d (deg) 




Count Rate (s 1 ) 


Flux f 


2785 


2002 May 13 


52407 


3.24 


8fl3 


55.7 


259.54098(8), -37.31437(6) 


99, 5 


0.0017(2) 


9(2) 


10131 


2009 Feb 19 


54881 


0.44 


0.'07 


32.0 


259.54098(1), -37.31419(1) 


81, 0.8 


0.0025(3) 


8(3) 


10766 


2009 May 15 


54966 


0.44 


0.'07 


33.3 


259.54096(1), -37.31439(1) 


82, 0.6 


0.0024(3) 


9(2) 


10767 


2009 Jul 28 


55040 


0.44 


0.'07 


34.2 


259.54088(1), -37.31432(1) 


66, 0.9 


0.0019(2) 


5(2) 


10768 


2009 Oct 23 


55127 


0.44 


0.'07 


34.1 


259.54093(1), -37.31432(1) 


73, 1 


0.0021(3) 


7(2) 



Frame time of the read out mode. 

The pointing offset from PSR J1718-3718 . 

Total live exposure time. 

Position of the X-ray counterpart of PSR .11718 — 3718 reported by the CIAD celldet 
nccrtaintics listed in this column were calculated using the source detection code, 
Total counts iV Brc and estimated background counts N^^^ in the 
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Fig. 1.— ATCA radio image of PSR J1718-3718 at 20 cm, show- 
ing the pulsar at its "on" phase and overlaid with the 2cr error el- 
lipse from the Parkes radio timing position (Manchester & Hobbs, 
2011 in preparation). The grey scale of the map is linear, ranging 
from —0.8 mjy beam -1 to +3.5 mjy beam -1 and the size of the 
restoring beam is shown at the lower left. The white spot marks 
the Chandra position, and its size is larger than the position un- 
certainty. 

rays that come from a point source of the same spec- 
trum as PSR J1718-3718(see ^T^. Then we em- 
ployed MARX to project the rays onto the detector where 
PSR J1718— 3718was located. For the above-mentioned 
images, we removed the effect of pixel randomization^ 
to improve their sharpness. The PSF broadening caused 
by the aspect reconstruction errors and ACIS pixeliza- 
tion were modeled by setting the DitherBlur param- 
eter to Of! 2 in MARX. We did not find any significant 
difference in the radial profile between the actual im- 
ages of PSR J1718— 3718 and the simulated images. We 
also aligned and merged all four PSR J1718— 3718 images 
from the 2009 observations to a single image, using 
the reproject_events and dmmerge tools in CIAD 4.2. 
Again, no significant difference was found between the 
point source's radial profile in the merged image and in 
the simulated image. In summary, we found no evidence 
of extended emission in the 2009 Chandra observations 
of PSR J1718-3718. 

In order to obtain a precise radio position of 
PSR J1718— 3718, we carried out a radio imaging cam- 



ct tool. Numbers in parentheses are ltr uncertainties in the last quoted digit. The 
ind are much smaller than the pointing uncertainty of Chandra . 
0.8-2.0 keV. 

fit of the spectra with JV R and kT fixed at their best-fit values. 

paign using the Australia Telescope Compact Array 
(ATCA) at 20 and 13 cm. A 12-hr observation was 
made on 2007 July 9 at center frequencies of 1384 MHz 
and 2368 MHz, with the 6C array configuration that 
gives the longest baseline of 6 km. During the obser- 
vation, pulsar gating was employed to record the pul- 
sar phase information. We performed the data reduc- 
tion using standard techniques in the MIRIAD pack- 
age. After the calibration process, a usable bandwidth 
of 104 MHz was left for each frequency band, which was 
split into 13 x 8 MHz channels. We then employed the 
task psrfix to de-disperse the pulsar signals and di- 
vided the data into 16 pulsar phase bins to form indi- 
vidual intensity maps separately. Our final maps have 
restoring beams of FWHMs 6'.'2 x 10"4 and 3'/2 x 5'.'7 
in the 20 cm and 13 cm wavebands, with corresponding 
rms noise of 0.5 mjy beam -1 and 0.6 mjy beam -1 , re- 
spectively. These are higher than the theoretical noise 
levels due to the sidelobes of a bright supernova rem- 
nant G349.7+0.2 at 8' to the South. The pulsar is 
clearly detected at the > 6<r level in the 20 cm wave- 
band (Figure [1]), but not at 13 cm. Finally, we employed 
the task imfit to determined the pulsar position, and 
found R.A.=17:18:09.84(5) and decl.=-37:18:52. 3(1.4) 
(J2000), with a flux density 3.5 ± 0.8 mjy in one phase 
bin (corresponding to a mean flux density of 0.22±0.05 
mjy). This flux density is consistent with that in the 
ATNF catalog for pulsars. This new radio position 
and the average Chandra X-ray position are only ~ 0'.'8 
apart, i.e. consistent with each other. 

2.2. Spectroscopy 

We extracted the spectrum of the pulsar from all five 
observations using the psextract script of CIAD 4.2. For 
the 2009 observations, we used a source region of radius 
3". A source region of radius 9'.'8 was used for the 2002 
observation because the source was offset from the center 
of the field of view, where the PSF is broader. The back- 
ground spectra were extracted from four circular regions 
of radius 10" centered around the pulsar for all the ob- 
servations. The resulting source and background spectra 
were then combined with RMF and ARF files generated 
using psextract and grouped with a minimum of 15 
counts per bin. We found that the spectra of the pulsar 
are soft, with very few counts above 2.0 keV. Thus for the 



' http: //exc .harvard. edu/ciao/why/acispixrand .html 



' http: //www. atnf . csiro . au/research/pulsar/psrcat/ 
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Energy (keV) 

Fig. 2. — Summed Chandra ACIS spectrum of all five observa- 
tions of PSR J1718— 3718 . The spectra are binned to contain a 
minimum of 20 counts per bin. The solid curve is the best-fit ab- 
sorbed blackbody model. The \ m the bottom plot is defined as 
the difference between the value of the spectral bin and the model 
prediction, divided by the uncertainty of the spectral bin. 

following spectral analysis we used only the 0.8-2.0 keV 
band. 

We fitted the five spectra separately with an absorbed 
blackbody model using XSPEC0 version 12.5.0. In each 
source spectrum, there are 66-99 total counts in 0.8-2.0 
keV (Table [T]), so the best- fit model parameters could 
not be well constrained when fitting one spectrum at 
a time. We fitted all five spectra jointly with a single 
blackbody model, and found a good fit with a reduced 
X 2 of 0.97 for 16 degrees of freedom. In a second joint 
fit, we allowed the normalization parameter to vary from 
observation to observation while fixing Nh and kT°° at 
their best-fit values, and found that the inferred 0.8-2.0 
keV absorbed fluxes were consistent with being constant 
(Table Q]) . This suggests that there are no statistically 
significant spectral or flux variations from observation to 
observation. Based on the 0.8-2 keV absorbed fluxes, we 
estimate a 3cr upper limit of 60% on any flux variations. 

Because the individual spectra have very few spectral 
bins after grouping, hence poor spectral resolution, they 
cannot constrain the model parameters well. In order to 
mitigate this problem, we summed the five spectra into 
a single spectrum. The resulting summed spectrum was 
grouped with a minimum of 20 photons per bin, and had 
18 spectral bins in 0.8-2.0 keV. The spectral resolution 
of the summed spectrum is much better than those of 
the individual spectra. 

We fitted the summed spectrum with a black- 
body model, neutron star atmosph ere model (NSA; 
iZavlin et all (l996t iPavlov et all 119951) and a power-law 
model, using the wabs model for interstellar absorption. 
Figure [2] shows a plot of the summed spectrum with the 
best-fit absorbed blackbody model. 

We found a best-fit blackbody temperature of 
186+ 1 geV, corresponding to a blackbody radius of 
1.8^0 5^4.5 km and a bolometric luminosity of 4+2 x 
lO 32 ^ 5 ergs -1 (assuming a fiducial distance d of 
4.5 kpc). In order to explore the confidence range of 
the redshifted temperature kT°° and radius for the 
blackbody model, we plotted their confidence contours in 
the left panel of Figure [3] This indicates the lowest pos- 

17 http : //heasarc . nasa . gov/docs/xanadu/xspec/ 



TABLE 2 

Spectral models for PSR J1718— 3718 and their best-fit 

PARAMETERS. 



Parameters Blackbody NSA a 



N H (10 22 cm" 2 ) 


1.3 ±0.2 


i 7 +o.i 

1 -'-0.2 


kT°° (eV) 




57l 7 2 


R°° (km) 




13 (fixed) 


Distance 13 (kpc) 




1 2+ 1 ' 4 
L - Z -0.7 


/ abs c (lO-^ergs- 1 cm" 2 ) 


0.078 ± 0.004 


0.077 ± 0.004 


/unab s d (lO-^ergs- 1 cm- 2 ) 


7 +0 - 4 
u - '-0.2 


0.80 ± 0.04 


L x c (10 32 erg s" 1 ) 


*_2 a 4.5 


n+3 
z -l 


xl(y) 


0.7(15) 


0.7(15) 



a The Hydrogen Atmosphere model for pulsar with B — 10 13 G 
and a pure hydrogen atmosphere. The values of neutron star sur- 
face temperature and radius R. — 10 km were redshifted for observers 
at infinite distance according to T°° = T(l - 2GM/Rc 2 ) 1/2 and 
R°° = R(l - 2GM/Rc 2 y 1/2 with M fixed at 1.4M . 
10 In the blackbody fit, a fiducial distance d of 4.5 kpc is used as a 
scaling factor for the best-fit parameters. In the NSA fit, distance is 
fitted. 

c Absorbed X-ray flux in 0.8-2.0 keV. 
d Unabsorbed X-ray flux in 0.8-2.0 keV. 

c Bolometric luminosity. For the NSA model, it is calculated based 
only on the pulsar's best-fit surface temperature, and docs not depend 
on distance. However, the best-fit distance in this fit is unreasonably 
small. If a more reasonable distance were assumed, the resulting best- 
fit temperature would likely become larger and thus imply a higher 
bolometric luminosity. 

sible kT°° of 140 eV, corresponding to «10km and 

L bol ~ 5xl0 33 ergs" 1 (higher than E = 1.6 xlO 33 ergs -1 
of the pulsar). 

The NSA model assumes that the X-ray emission of 
the pulsar comes from its entire surface. The best-fit 
local surface temperature kT is only 75~tio e ^ (kT°° — 
57t 7 2 eV as seen from Earth), with a best-fit distance 
of 1.2+0 7 kpc ( mucn smaller than the estimated range 
of 4.5-10 kpc) assuming that the neutron star has mass 

1.4 Mq and local radius 10 km. The best-fit parameters 
of the blackbody and NSA models are listed in Table 
[21 The highly magnetized NSA model assumes a B field 
of 10 13 G, less than the inferred B of 7.4 x 10 13 G for 
PSR J1718— 3718 . We allowed the normalization param- 
eter, which corresponds to 1/d 2 , where d is the distance 
of the neutron star, to vary when fitting the spectrum. 
We plot the confidence contours of the redshifted effec- 
tive surface temperature kT°° and distance in Figure [31 
right panel. Assuming the pulsar is at a distance between 

4.5 kpc and 10 kpc, it should have a surface temperature 
between 75 eV and 97 eV and a bolometric luminosity 
> 10 33 ergs" 1 . 

The best-fit absorbed power-law model has an unrea- 
sonably large photon index (> 8); therefore, we con- 
sider it no further. We also tried to fit the 0.8-10.0 keV 
summed spect rum with a reso nant cyclotron scattering 
model (RCS) (|Rea et all 120081 ) . but did not find a good 
fit (best reduced \ 2 — 2.4 for 15 degrees of freedom); 
this is likely due to the lack of hard photon events in the 
source spectrum. 

2.3. Variability and Pulse Profile 

We adjusted the time stamps of the source events from 
all five observations to the solar system barycenter time 
using the axbary tool in CIA0. We binned the photon 
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Fig. 3. — Left panel: 68%-, 90%- and 99%-confidence contours of R^ (assuming a fiducial distance of 4.5 kpc) and kT°° from blackbody 
fitting. Right panel: confidence contours of pulsar distance and fcT°° from fitting with the NSA model with B = 10 13 G, assuming the 
neutron star radius R = 10 km (R°° = 13 km). The dashed lines in both panels are the contours of constant bolometric luminosity in units 
of ergs -1 . The shaded area in the right panel marks the region for which the pulsar's distance is 4.5—10 kpc. 
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Fig. 4.— Pulse profiles for PSR J1718-3718 in the 0.8-2.0 keV 
energy range. Also given in the plot are the reduced x 2 values for 
7 degrees of freedom from fitting the profile with a constant. The 
peak of the pulsar's radio pulse aligns with zero phase in this plot. 

events of energy between 0.8 and 2.0 keV from the five ob- 
servations evenly in time with 3.4 hr per bin. The result- 
ing count rates were consistent with being constant with 
a 3cr upper limit of 48% on variations, therefore showing 
no evidence of significantly variability on timescales of 3- 
9 hr. In order to look for variability on shorter timescales, 
we measured the intervals between the arrival times of 
every two photons in each observation. We found that 
they are consistent with the exponential distribution ex- 
pected from a constant count rate, and therefore, show 
no evidence of flux variations. 

Unlike the 2002 normal timed exposure mode obser- 
vation that has time resolution 3.24s, the later four 1/8 
subarray observations have 0.44 s time resolution (Table 
[1]) and therefore could be used to search for pulsations 
from this 3.3-s pulsar. We folded the events with ener- 
gies between 0.8 and 2.0 keV from the 2009 observations 
into 8 phase bins based on a timing ephemeris obtained 
using the Parkes telescope (Manchester & Hobbs, 2011 
in preparation). The resulting pulse profile is shown in 
Figure HJ Significant pulsations are detected. We found 



TABLE 3 

Significance of the X-ray pulsations of 
PSR J1718-3718. 



Observation 


X? 


p x 2 

null 


H value 


pH 
^iiull 


Summed 


49.0 


2 x 10~ 8 


44.5 


2 x 10~ 8 


2009 Feb 19 


6.2 


0.5 


3.0 


0.3 


2009 May 15 


10.1 


0.2 


8.1 


0.04 


2009 Jul 28 


36.6 


6 x 10 -6 


24.4 


6 x 10~ 5 


2009 Oct 23 


16.2 


0.02 


11.8 


0.009 



the H test (|de Jagerlll994D value of the summed profile 
is 44.5 and the best-fit reduced \ 2 is 7.0 for 7 degrees 
of freedom. Both correspond to null-hypothesis possibil- 
ities of ~ 2 x 10~ 8 , clearly excluding the null hypoth- 
esis . We also measured an area pulse d fraction (the 
fracti onal counts above the minimum; Gonzale z et al.l 
12010( 1 of 52%±13% in the 0.8-2.0 keV band and a max- 
min pulsed fraction (N max - N min )/(N max + N min ) of 
60%±13%. 

In Figure SI we also plot the folded pulse profiles from 
individual observations. Not all of them were signifi- 
cantly pulsed. For instance, when fitted with a constant, 
the pulse profile of the 2009 February 19 observation 
gives a best-fit reduced x 2 01 0-9 f° r 7 degrees of free- 
dom and a H value of 3, which do not exclude the null 
hypothesis (TableEJ. Through numerical simulations, we 
have verified that with only 81 counts and assuming Pois- 
son noise, it is possible for a source having 52% pulsed 
fraction to produce a pulse profile of such low signifi- 
cance. We simulated 10000 pulse profiles with a source 
with area pulsed fraction of 52%, and found that 209 of 
them show lower pulse significance than in the February 
19 observation. Thus, even if the pulsar's profile did not 
change between 2009 February 19 and July 28, there is 
~2% chance of observing a pulse profile similar to the 
February 19 one. Taking the number of trials into ac- 
count, the low pulse significance of this observation does 
not provide strong evidence for a change in the pulsed 
fraction. 

We fitted the X-ray pulse profile with a sinusoidal func- 
tion to find its peak phase and compared it with that of 
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the radio pulse after correcting for the effect of disper- 
sion due to the interstellar medium. The radio pulse 
leads X-ray pulse by 0.01 ±0.03 in phase. Thus, they are 
consistent with being aligned. 

3. DISCUSSION 

We have reported on four new Chandra observations, 
plus a reanalysis of one archival observation, of the 
young, high-B RPP J1718-3718. We found no evidence 
of magnetar-like flux variability in PSR J1718— 3718 from 
our Chandra observations, and set a 3d upper limit on 
any flux variability of 60% in the 0.8-2 keV band. How- 
ever, the possibility that a magnetar-like outburst, such 
as that observed from PSR J1846-0258 in 2006 which 
lasted for only few weeks, happened in the span of our 
observations could not be ruled out. 

PSR J1718— 3718 's X-ray spectrum is soft and ther- 
mal, and is well fit by a blackbody model. Fitting its 
summed spectrum with a blackbody model, we found a 
high blackbody temperature of 186l 1 geV (Tabled]), and 

a corresponding best-fit blackbody radius of l.Sij^di.s 
km. Such an emission radius is consistent with radiation 
from hot spots. However, it is not consistent with polar 
caps heated by return currents because of the unusually 
high X-ray efficiency (L^/E = 0.3d| 5 ). By contrast, 
models for polar-cap heating predict that no more than 
~ 10~ 3 of the spin-dow n luminosity should be con verted 
to thermal radiation (jHarding fc Muslimov|[200~ll ). Note 
that if the distance is larger than 4.5 kpc, this conclusion 
is only strengthened. Indeed at 10 kpc, > E. Based 
on the confidence contours of kT°° and (left panel 
of Figure [3]), we cannot completely exclude a blackbody 
fit of kT°° = 140 eV and R^ h = 10d 4 . 5 km. However, 
the measured 52%±13% area pulsed fraction suggests 
that the surface temperature of the pulsar cannot be uni- 
form. Given the pulsar's spin-down age of 34kyr, a sur- 
face temperature of 140 eV is still higher t han what one 
would expect (60-90 eV. iPage et all 12001 ) from a mini- 
mum cooling model for the neutron star surface without 
considering the effects of the magnetic field. Interest- 
ingly, the 186 eV best-fit blackbody temperature is sim- 
ilar to those found for the high-B PSRs J1119-6127, 
J1734-3333, J1819-1458 (see TableHfor details and ref- 
erences) and the transient AXP XTE J18 10-197when it 
was i n quiescence between 1980 and 1993 (Go tthelf et al.l 
l200l kT°° = 180 ± 10 eV). 

On the other hand, fitting the spectrum with a NSA 
model leads to a best-estimated surface temperature of 
75-97 eV (assuming a neutron star mass of 1.4M , a local 
radius of 10km, and a pulsar distance of 4.5-10 kpc), 
consistent with standard cooling. We note, however, that 
the magnetic field strength assumed in the NSA model is 
10 13 G, almost one order of magnitude smaller than the 
spin-down-inferred value. Therefore, the results of the 
NSA model fit should be taken with caution. 

Attempting to explain the X-ray thermal emission 
observed from magn e tars, X INSs and so me h igh- ,8 
pulsars. I Arras et all (|2004l ). iPons et al.l (|2007l ) and 



lAguilera et al.l ( 2008 ) constructed neutron-star cooling 
models in which pulsars with magnetic fields higher 
than 10 13 G are significantly heated by field decay. 
The key evidence to support this theory is an in- 
triguing possible correlation found between the pul- 
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Fig. 5. — Blackbody temperatures versus characteristic ages 
(P/(P(n — 1)), where n is the measured braking index, assum- 
ing n = 3 if it is not measured) of high-B pulsars (filled circles), 
normal pulsars (open triangles), and XINSs (filled squares). Ref- 
erences are listed in Tableful 

sar's blackbody temperature T and spin-down magnetic 
field B (T oc B 1 / 2 ) based on a sample of magne- 
tars, XINSs and some RPPs (jPons et al.l 120071 ). How- 
ever, their analysis did not consider high-_B RPPs. 
Searches for evidence of magnetic-field-decay heating 
have been conduct ed on several high-i? pulsars, e.g., 
PSRs J1814-1744 dPivovaroff et al.ll2000h. J1847-0130 



dMcLaughlin et al.l I2003D. B0154+61 dGonzalez et all 
I2004T). J1119-6127 (iGonzalez et all 120051) . J1718-3718 
dKaspi fc McLaughlin! I2005D. B1916+14 (|Zhu et all 



|2009j) and J1734-3333 (|Qlausen et aHl2010h . and X-ray 
counterparts were found in some cases. Their spectra, 
however, have not yet been sufficiently well constrained 
to prove the existence of significant magnetic-field-decay 
heating. For a recent review on high-magnetic-field pul- 
sars, see INg fc Kaspil (120101) . 

With our new data and spectral analysis of 
PSR J1718— 3718, we are unable to confirm that it is 
heated by magnetic field decay. This is mainly because 
the non-magnetized neutron star cooling models predict 
a large range of surface temperatures for a given pulsar 
age. However, we can compare the surface temperatures 
of several high-i? RPPs with those of normal RPPs, pro- 
vided that the temperatures were measured using the 
same spectral model In Figure [5j we plot blackbody 
kT°° versus age for a collection of pulsars including some 
high-_B pulsars (see Tabled]). From this plot, one can see 
that the blackbody temperatures of the high-S pulsars 
appear to be in general higher than those of the normal 
pulsars. 

We also looked for the same T-B correlation showed 
bv IPons et all (|2007f ) in a kT°° versus B plot, but the 
temperatures of the pulsars are too scattered to discern 
a trend. This could be because our sample has a small 
range of B but a large range of ages. 

Note that three pulsars listed in Table 4 (PSRs 
B1929+10, B0355+54 and J1357-6449) are not plot- 
ted in Figure |3J This is because they all exhibit a large 
kT°° with very small blackbody radius, R^ < 1 km, 
consistent with return-current heating. For the other 
RPPs, we cannot rule out the possibility that their 
blackbody temperatures are also higher because of re- 
turn currents, but these sources nevertheless provide in- 
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TABLE 4 

Surface temperatures measured for high-B pulsars, normal pulsars, and XINSs. 



PSR 


T c (kyr) 


B (G) 


fcT-(eV) 


i?-(km)[D(kpc)] 




B0950+08 


18000 


O A \s 

Z.4 X 


mil 

1U 


<41.0 


10[0.3] 


<0.06 


B1929+10 c 


3100 


5.2 x 


i nil 
10 


300 + ?n 


0.033+S?,?,!? [0.41 

— 0.005 1 J 


3 x 10~ 4 


J0538+2817 


40 


7.3 x 


10 11 


181 ± 3 


2.23 ± 0.01[1.5] 


0.01 




564 


8.4 x 


10 11 


200+ 96 


n i2 +0 - 16 n ol 

u ' lz — 0.07 l 1 '"! 


7 x 10 — 5 




^yuu 


9.6 x 


10 11 




i nrn ^1 

lU[U.OJ 


in 

<,U. 1U 


J_j J. WUtJ — 


535 


1.1 X 


i nl2 
10 


68 i 3 


l^-o_0.71 u ' ,5 J 


0.01 


J0633+1746 


342 


1.6 x 


10 


41.4 ± 0.1 


9 ± 1[0.2] 


9 x 10~ 4 


J1811-1925 


23 


1.7 x 


i nl2 
10 


<150 


10[5.0] 


<0.001 


11 740+1 oon 


114 


1.8 x 


i nl2 
10 


70+ 


7 Of! 41 


7 x 10~ 4 


R1 823—1 3 


21 


2.8 x 


10 12 


97+4 


6.3[4.0] 


2 x 10~ 4 


R1 7f\f{ — AA 
J_J J. i UvJ 'irfc 


18 


3.1 x 


10 12 


143 ± 14 


q « 4. n QT9 ^1 


2 x 10 — 4 


DflOOO A E 

rSUooo — 40 




3.4 x 


10 12 


no _l_ o 


c l+0.4r n ol 


A \s 1 fl — 5 


B1046— 58 


20 


3.5 x 


10 12 


<95.0 


10[2.7] 


<5 x 10 — 4 


TnoriK i a a a n 
J 0205+6449 


2.4 


3.6 x 


10 12 


112 i 9 


11[3.2] 


9 x 10 


T3nc;Q 1 i 01 
L5Uool+zl 


U.yo 


3.8 x 


10 12 


<1 i 1 


16[1.7] 


r \y 1 n — 5 
<o X 1U 


B0656+14 


111 


4.7 x 


10 12 


56.0 ± 0.9 


2lt 3 [0.3] 


0.01 


J1357-6429 c 


7.3 


7.8 x 


10 12 




1.0[4.1] 


3 x 10~ 5 


B2334+61 


41 


9.9 x 


10 12 


109 ± 35 


1.7[3.1] 


8 x 10~ 4 


J1856-3754 b 


3800 


1.5 x 


10 13 


63.5 ±0.2 


6.2±0.1[0.2] 


24 


B1916+14 


88 


1.6 x 


10 13 


130tr 


0.8±0.1[2.1] 


0.005 


J2143+0654 b 


3700 


2.0 x 


10 13 


104 ±4 


3.1[0.4] 


76 


B0154+61 


197 


2.1 x 


10 13 


<73.0 


10[1.7] 


<0.6 


J0720-3125 b 


1900 


2.5 x 


10 13 


90 ±4 


6.4[0.4] 


73 


J0806-4123 b 


3300 


2.5 x 


10 13 


87 ± 11 


1.3[0.2] 


8 


J Uo't ( - — ^olO 


7Qfl 


9 7 v 








^ c; o 
\ Oo 


J1846-0257 


442 


2.7 x 


10 13 


<120 


10[5.2] 


<38 


J1308+2127 b 


1500 


3.4 x 


10 13 


100 ±2 


4.1[0.5] 


54 


J1119-6127 


1.8 


4.1 x 


10 13 


210 ± 10 


2.7 ± 0.7[8.4] 


8 x 10~ 4 


J0420-5022 b 


109 


4.2 x 


10 13 


45 ±3 


3.3[0.3] 


0.01 


J1846-0258 


0.88 


4.9 x 


10 13 


<250 


2.7[6.0] 


<5 x 10" 4 


J1819-1458 


117 


5.0 x 


10 13 


120 ± 20 


2.1 ± 0.4[3.6] 


0.4 


J1734-3333 


8.1 


5.2 x 


10 13 


2501"° 


1± 3 [6.1] 


0.01 


J1814-1744 


85 


5.5 x 


10 13 








J1718-3718 


34 


7.4 x 


10 13 


i89±y 


1-8±J:J[4.5] 


0.3 


J1847-0130 


83 


9.4 x 


10 13 








The ratio of the pi 


llsar's boloi 


nctric h 




ity (tgf, s 1.28 


X 10 35 (R~ ) 2 (fcT) 4 erg a 


- 1 ) to spin dow, 
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McGowan et al. (20071 
Becker et al. (20041 
De Luca et al. (20051 
De Luca et al. (2005) 



Kaspi et al. (2006 



Z. Misanovic ct al. (2011, in prep.) 

Pavlov et al. (2008) 

Gotthelf et al. (2002) 

Manzali et al. (2007) 
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Weisskopf et al. (20041 
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Rea et al. (2009a) 

Olausen et al. (2010) 
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XINSs. 

C These pulsars have very small blackbody radius < 1km and L^/E < 10~ 3 . Their thermal radiation is likely coining from hot spots 

by return current heating. Therefore they arc not included in Figure fol 



teresting upper limits on any non-return current ther- 
mal emission. On the other hand, return-current heat- 
ing is unlikely to be present from the high-B PSRs 
and XINSs. This is because most of them, including 
PSR J1718-3718, (B > 10 13 G, see Table g]) have an 
X-ray efficiency L^/E > 1CT 3 . Such a high X-ray ef- 
ficie ncy is clearly inconsistent w ith return-current heat- 
ing (jHarding fc Muslimovl 120011 ). Also note that the X- 
ray spectra of some pulsars, such as PSRs B1055— 52, 
B0656+14 and J0633-I-1746, show evidence of thermal 
emission from both a hot spot and a much cooler neu- 
tron star surface; in these cases, only the kT^ of the cool 
surface was used. 

In summary, our Chandra observations of 
PSR J1718— 3718 have revealed, for the first time, 
X-ray pulsations at the pulse period, as well as a 
thermal spectrum of blackbody temperature somewhat 
higher than for other RPPs having the same age. We 
have found a high bolometric to spin-down luminosity 
ratio, ~0.3 for a distance of 4.5 kpc, and higher for 
more realistic, larger distances. Although we cannot 
rule out standard passive cooling, as a model fit with an 



NSA model yields a lower surface temperature, we have 
considered the possibility that PSR J1718— 3718 exhibits 
enhanced thermal emission due to magnetic-field decay, 
as predicted by m o dels of mag neto-thermal evolution 
(I Arras et all 120041: IPons et all 120071: lAguilera et all 
I2008D . We have compiled such measurements for the 
other high- and low-i? RPPs, and find a hint that those 
of higher B are generally hotter than Iow-jB pulsars of 
the same age. However deeper observations of high- and 
low-B pulsars are required to confirm this possibility. 
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